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Slice-stretching effects are discussed as they arise at the event horizon when geodesically slicing 
the extended Schwarzschild black-hole spacetime while using singularity excision. In particular, 
for Novikov and isotropic spatial coordinates the outward movement of the event horizon ("slice 
sucking") and the unbounded growth there of the radial metric component ("slice wrapping") are 
analyzed. For the overall slice stretching, very similar late time behavior is found when comparing 
with maximal slicing. Thus, the intuitive argument that attributes slice stretching to singularity 
avoidance is incorrect. 
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I. INTRODUCTION 

When evolving a spacetime containing a physical sin- 
gularity without making use of a shift, the foliation 
usually is of pathological nature as so-called "slice- 
stretching" effects develop pj . Here "slice sucking" arises 
in the form of outward-drifting coordinate locations as 
the corresponding observers arc falling towards the sin- 
gularity. This infall takes place in a differential manner, 
leading to large proper distances in between neighboring 
observers with "slice wrapping" showing up in the form 
of large gradients in the radial metric function. 

In the past, slice-stretching effects have been often at- 
tributed to the singularity-avoiding behavior of a folia- 
tion by the following intuitive argument: For such slices 
the lapse collapses to zero in the strongly curved "inte- 
rior" region and the evolution of the metric essentially 
freezes there, while it marches ahead further "out" in or- 
der to evolve a large fraction of the spacetime. Hence, 
singularity-avoiding slices "wrap up around the singu- 
larity" 13 which causes "large amounts of shear in the 
coordinate grid" Q. 

To study whether the singularity- avoiding property of 
a foliation plays a role for the overall slice stretching, I 
will in this paper for the extended Schwarzschild space- 
time compare slice-stretching effects for geodesic and 
maximal slicing. 

Characterized by unit lapse and vanishing shift, 
geodesic slicing represents the "simplest" gauge choice 
which, however, does not avoid physical singularities. 
In particular, starting with the time-symmetric con- 
formally flat Einstein-Rosen bridge 3 of the extended 
Schwarzschild spacetime, an observer at the throat is ini- 
tially at the event horizon teh = 2M and falls freely into 
the singularity. Hence, if the singularity is part of the 
grid, a numerical simulation faces infinite curvature and 
has to crash after evolving for the free-fall time given by 
ttM. This fact has been used frequently in numerical 
relativity when performing the so-called "crash test" as 
e.g. in kj or [BJ. In addition, when not excising the sin- 



gularity, the analytically known non-trivial evolution of 
the 3-metric can be used up to the time ttM for testing 
numerical codes and for studying the numerical behavior 
of different formulations of the Einstein evolution equa- 
tions, see e.g. Q and 0- 

In the following, however, I am only interested in the 
portion of the geodesic slices lying in the exterior parts 
of the spacetime. From the viewpoint of numerical rela- 
tivity, I make use of a code capable of excising the inte- 
rior region in between the "left-hand" and "right-hand" 
event horizon. For Novikov spatial coordinates this idea 
is shown in Fig. I of Q. In the same reference it has 
also been demonstrated numerically that geodesic slicing 
together with singularity excision can be used to evolve 
a single black hole for considerably more than ttM. 

When studying the late time behavior of the geodesic 
slices, I will for simplicity concentrate on the event hori- 
zon acting as a "marker" for slice-stretching effects. In 
particular, the location of the event horizon in terms of 
Novikov and isotropic spatial coordinates will be deter- 
mined as a function of time together with the behavior 
of the radial component of the 3-mctric there. 

These results will then be compared both analytically 
and numerically with statements obtained for maximal 
slicing in [2, flOt Hll Il2l ] . Being motivated geometrically 
and characterized by the condition that the trace of the 
extrinsic curvature vanishes at all times fl3T |. maximal 
slices avoid the Schwarzschild singularity by approach- 
ing the limiting slice r = 3M/2 asymptotically [14(. The 
main result of this comparison is that the overall slice 
stretching at the event horizon for both geodesic and 
maximal slicing in leading order is found to be propor- 
tional to time. Hence the intuitive argument that at- 
tributes slice stretching to singularity avoidance turns out 
to be incorrect. 

The paper is organized as follows: In Sec. [H] slice- 
stretching effects at the event horizon are discussed, 
studying in Subsec. Ill Al the cycloidal motion of ra- 
dial geodesies and focusing in Subsecs. Ill Bl and III CI on 
Novikov and isotropic spatial coordinates, respectively. I 
conclude in Sec. II I II 
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II. SLICE STRETCHING AT THE EVENT 
HORIZON 

A. Cycloidal motion of radial geodesies 

Geodesic slicing is characterized by unit lapse and van- 
ishing shift which define Gaussian normal coordinates. 
Those are comoving in the sense that radially freely 
falling observers are at rest and the time coordinate mea- 
sures proper time. One important property of Gaussian 
normal coordinates is that the geodesies defining the co- 
ordinates remain orthogonal to all constant time hyper- 
surfaces, and transformations between different spatial 
coordinates are hence time-independent. 

In particular, referring to Schwarzschild coordinates, a 
radial geodesic starting at the singularity at r = per- 
forms a cycloidal motion out to some maximal radius 
f and back as pointed out in more detail in §25.5 and 
§31.3 of [13 • in terms of r > 2M, the geometry of the 
Schwarzschild spacetime is described by the line element 



ds 2 = -dr 2 
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Here by integrating the geodesic equation one can see 
that r = r(r, r) is given implicitly by 
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Furthermore, following it turns out that by implicit 
differentiation 
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is found (but note a missing root in the formula given in 
that reference). 

Concentrating for a discussion of slice stretching on 
the event horizon, r = teh = 2M, in leading order (to 
be denoted by ~) from one can infer 
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Interested in the behavior of the radial component Gff 
of the line element Q there, when using © at late times 
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is found. 

In a next step it is then desirable to combine the 
coordinate-dependent effects of slice sucking and slice 
wrapping, Q and iJSJ, into one coordinate- independent 
quantity measuring the overall slice stretching. Here I 
want to point out that for maximal slicing in |l2j the 



overall slice stretching has been characterized in a global 
sense by integrals of metric quantities such as the square 
root of the radial metric component. Those integrations 
in space have been carried out over the throat and, when 
integrating e.g. from the left- to the right-hand event 
horizon, for "numerically favorable" boundary conditions 
a divergence being in leading order proportional to time 
has been found, see ^| f° r details. Note, however, that 
one can not carry out such an integration over the throat 
for slicings like geodesic slicing which do not avoid but 
hit the singularity. 

Therefore I now want to characterize the overall slice 
stretching in a local sense only. For this task I introduce 
the corresponding indefinite integrals and evaluate them 
at the right-hand event horizon setting r = teh = 2M 
and using previously derived results. Integrating the 
square root of the radial metric component in this sense, 
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for geodesic slicing one obtains a divergence which is in 
leading order linear in time. 

For maximal slicing one can write the radial met- 
ric component as G xx = x A ^ 12 {x) / r A . This formula 
was found in [10| for evolutions using isotropic grid 
coordinates x where the conformal factor is given by 
fy(x) = 1 + M/2x. In terms of these coordinates the lo- 
cation of the right-hand event horizon is described by 
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as pointed out in |l2j| . and hence with 
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also for maximal slicing a behavior of order C(r) is ob- 
tained [It| . 

The coordinate-independent observations ijHJ and (jHJ 
clearly are in contradiction to the argument that slice 
stretching is due to singularity avoidance. In order to 
make further statements regarding the "splitting" of the 
overall slice stretching into slice sucking and slice wrap- 
ping, I will discuss in the next subsection geodesic slicing 
in terms of (more familiar) Novikov coordinates. 



B. Novikov coordinates 

The Novikov spatial coordinate R* is related time- 
independently to f by 



2M{R* 2 + 1) 



(9) 
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where the isometry R* < — > —R* is mapping the two 
"universes" of the extended Schwarzschild spacetime into 
each other. Considering Q), one can then observe that 
in terms of R* slice stretching takes place in a symmetric 
manner as the locations of the left- and right-hand event 
horizon at late times are given by 



R 



EH 
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(ttM)3 

For the radial component of the metric, 
G r * r *(t,R*) = 16M 2 (R* 2 + 1) 



(10) 



(11) 



using © in leading order at both left- and right-hand 
event horizon an unbounded growth like 



Gr*r* 



9(ttM) 4 < 



(12) 



is found. 

I want to emphasize here that slice sucking and slice 
wrapping at the event horizon - for Novikov spatial co- 
ordinates - tak es p lace in order ©(r 1 / 3 ) and C(r 4 / 3 ), re- 
spectively. In the same late time behavior - but for 
isotropic grid coordinates - has been found for maximal 
slicing (at the right-hand event horizon in the context of 
even or "zgp" boundary conditions). 

In order to compare the slice-stretching effects for the 
same choice of spatial coordinates, in the next subsection 
I will geodesically slice puncture data. 



C. Isotropic coordinates 

Isotropic spatial coordinates can be introduced again 
in a time-independent way by 



(13) 



making use of the conformal factor ty(x) = l + M/2x. 
Since the isometry x < — ► M 2 /Ax is present, with Q and 
(I13|l one may readily verify 



M 2 



271/3 



4.c 



-T3 



(14) 



EH 



and observe that in leading order the location of the 
right-hand event horizon grows in order 0(t 2 / 3 ) whereas 
the left-hand event horizon approaches the puncture like 
0(r~ 2//3 ). Furthermore, it turns out that the radial com- 
ponent of the 3-metric - often rescaled by \& 4 (x) to focus 
on the dynamical features in the metric rather than on 
the static singularity at x = - is given by 
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When using (© at the left- and right-hand event horizon 
for Qxx an unbounded growth of order C(r 2 / 3 ) is found, 
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whereas when including the conformal factor for G xx a 
behavior of order 0(t 10 / 3 ) and 0(t 2 / 3 ) is obtained, 
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The geodesic slicingof puncture data has been investi- 
gated numerically in Here I show the corresponding 
spacetime diagram in Fig. ^ 
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FIG. 1: The geodesically sliced Schwarzschild spacetime in 
isotropic coordinates is shown, denoting the location of the 
left- and right-hand event horizon by downward- and upward- 
pointing triangles, respectively. Note that in numerical prac- 
tice one would usually excise from the puncture at x = up 
to a "ghost zone" to the left of the right-hand event horizon. 

Furthermore, by making use of the regularized spher- 
ically symmetric code described in 0], I have evolved 
Schwarzschild puncture data numerically with geodesic 
slicing and, for comparison, maximal slicing. For sim- 
ulations using 30,000 grid points and a resolution of 
Ax = 0.001M the location of the right-hand event hori- 
zon and the value of the radial metric component there, 
x~t M and G X x\„+ i are shown as a function of time in 
Fig. El 

One should observe in this figure that for geodesic slic- 
ing the numerically found slice sucking and slice wrap- 
ping is up to r = 10M in excellent agreement with the an- 
alytically predicted results when excising from the punc- 
ture up to 0.98 times the value of x~^ H . Due to growing 
errors at the excision boundary, however, the run fails 
shortly afterwards. 

Foliations of Schwarzschild puncture data using max- 
imal slicing, its lapse arising from the elliptic equation 
Aa = Rot, have been studied both analytically and nu- 
merically in H E3, [H Il2^ . In particular, for runs de- 
manding symmetry with respect to the throat or for 
the puncture evolution imposing a vanishing gradient of 
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FIG. 2: Slice sucking and slice wrapping are shown as tak- 
ing place for isotropic coordinates at the event horizon. Here 
the numerically obtained data for geodesic slicing (shown as 
stars) is in excellent agreement with the analytically predicted 
line, with both x%„ and G xx \ + being of order C(r 2//3 ). 

an ' x eh 

For maximal slicing these slice-stretching effects are of order 
0(t 1//3 ) and 0(t 4,/3 ), obtaining numerically very similar re- 
sults when demanding even or "zgp" boundary conditions 
(denoted by upward- or downward-pointing triangles, respec- 
tively) . 



the lapse at the puncture, i.e. for even or "zgp" bound- 
ary conditions, slice sucking and slice wrapping at the 



right-hand event horizon have been found to be of order 
©(t 1 / 3 ) and C(r 4 ' 3 ), respectively. Numerical results are 
shown here in Fig. 



III. CONCLUSIONS 

I have studied the slice-stretching effects which 
are present when geodesically slicing the extended 
Schwarzschild spacetime while making use of singular- 
ity excision. The analysis has been carried out at the 
event horizon in terms of Novikov and isotropic spatial 
coordinates. 

Independent of the coordinate choice, the overall slice 
stretching has been found to be proportional to time and 
hence to be comparable to the one arising for maximal 
slicing. Its splitting into slice sucking and slice wrapping, 
however, for the same choice of spatial coordinates turned 
out to be different. 

The intuitive argument that attributes slice stretching 
to singularity avoidance has been found to be incorrect. 
Instead, for evolutions with vanishing shift, slice sucking 
and slice wrapping are caused by the differential infall of 
Eulcrian observers. 

Taking the maximally sliced Schwarzschild spacetime 
as an example, it is furthermore of interest to study ana- 
lytically whether these effects can be avoided by making 
use of a geometrically motivated shift. I will report on 
work in this direction in a further paper [l9l | . 
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